Search for WIMPs in liquid argon§ 

c. amslerEI 

Physik-Institut der Universitat Zurich 
CH-8057 Zurich, Switzerland 

Abstract: 

Our group from the University of Zurich is performing R&D work towards the design of a large 
liquid argon detector to detect Weakly Interacting Massive Particles (WIMPs). This project is 
developed within the DARWIN Collaboration funded by ASPERA to prepare a proposal for the 
next generation of WIMP searches using noble liquids. We are performing R&D to detect the 
VUV light from recoiling argon nuclei. Results obtained with one ton of liquid argon (ArDM 
prototype) and prospects using a monoenergetic neutron source are discussed. 

1 Introduction 

From astronomical observations it is known that most of the gravitational mass in the universe 
does not couple to electromagnetic radiation and is therefore not directly observable with tele- 
scopes. Recent measurements lead to the conclusion that about 72% of the mass density is due 
to dark energy, 23% to dark matter (DM), and only 5% to baryons The contribution from 
neutrinos is less than 1%. Baryons are mostly located in the intergalactic gases with only about 
10% contributing to star masses. Roughly 15% of the DM may be due to MACHOS or cold 
molecular clouds. Although exotic mecanisms (such as modifications of Newton's gravitational 
law [2]) that do not require DM have been postulated, the natural assumption is that DM indeed 
exists as some form of elementary particles, axions or Weakly Interacting Massive Particles 
(WIMPs), or possibly as primordial black holes (having formed before nucleosynthesis). The 
direct laboratory observation of dark matter (DM) is one of the most pressing issues in Particle 
Physics, which can be addressed at the LHC and with non-accelerator experiments. 

Because DM has survived since the birth of the universe it has to be stable and only weakly 
interacting. As mentioned already, the most popular candidates for DM are axions and WIMPs. 
Heavy neutrinos would qualify as WIMPs, although there is no obvious reasons as to why 
they should be stable. The lightest superparticle in SUSY models conserving i?-parity is the 
most popular WIMP candidate, the spin 1/2 neutralino Xi with mass in the 10 GeV to 10 TeV 
range. i?-parity conservation ensures that the Xi is stable. Also, the Xi cannot transform into 
other SUSY particles when interacting with matter, due to its low mass. At the LHC the Xi 
will therefore manifest itself by a large missing energy. On the other hand, the can scatter 
e.g. on constituent quarks in nucleons or nuclei (Fig. [T] left), leading to nuclear recoils in 
the range of 1 - 100 keV. Non- accelerator laboratory searches are all based on the detection 
of such nuclear recoils. The scattering cross section on nucleons is tiny, in the range 10^^ to 
10^^ pb, comparable to that for neutrino interactions. The differential cross section decreases 
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exponentially with recoil energy (see Fig. |3] below) which makes WIMP detection difficult 
due to the low energy background. Hence massive detectors with low detection thresholds are 
needed, among them cryogenic ones using liquid argon (LAr) or liquid xenon (LXe). 
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Figure 1 : Left: Feynman graphs of x = Xi interactions with quarks in the nucleon. Right: 
WIMP-nucleon cross section as a function of WIMP mass showing the most stringent upper 
limit (90% CL)from XENON-100 (adapted from ref. [3]). Theoretical predictions from super- 
symmetry 14] are also shown. 

WIMPs can annihilate in the galactic halo, leading to the production of e.g. positrons, 
antiprotons, neutrinos or photons. For instance, an excess of 7-rays consistent with 80 GeV 
WIMPS has been reported in a re-analysis of EGRET data [5|. Satellite experiments such 
as GLAST, AMS2, PAMELA, Cerenkov detectors such as MAGIC and VERITAS, or large 
neutrino telescopes such as ICE CUBE, might reveal signals from WIMP annihilations |[T1. An 
annual modulation {63a effect) was reported by the DAMA experiment with 100 kg of Nal 
detectors [6 |, and more recently confirmed by the same group (8.2cr effect with 250 kg of Nal). 
The signal was, however, not observed by other experiments. In particular, the spin-dependent 
signal was recently ruled out by COUPP [7 |. 

Fig. [H (right) shows upper limits for WIMPs from recent experiments. Since one does not 
know precisely the cross section, one usually shows the limits in a two dimensional plot cross 
section vs. mass. For low masses the sensitivity decreases, due to the low recoil energy and 
the detection threshold. For high masses the loss of sensitivity is due to the diminishing WIMP 
flux. The best upper limits are obtained for a WIMP mass comparable to the atomic number A 
of the target. The spin independent rate provides the current best upper limits for WIMPs and 
is also the relevant interaction for our developments with liquid argon since ^°Ar is a spin zero 
nucleus. The XENON-100 [3] and CDMS [8] experiments have produced the best upper limits 
so far, the former even reaching 7 x 10^^ pb for masses around 50 GeV. The current upper 
limits now reach into the relevant region of theoretical predictions based on supersymmetry [4J. 

Most WIMP searches can be classified in two classes, (1) those using pure semiconductors 
and, (2) those using noble liquids such as LXe or LAr. In semiconductors operating at very 
low temperatures one measures the phonons and/or the ionization generated by recoiling nu- 
clei. Representatives of this class are CDMS, CRESST, EDELWEISS. In noble liquids (ArDM, 
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DEAP, WARP, XENON, ZEPLIN) one mesures the UV Hght and possibly also the electric 
charge produced by the ionizing recoils. Sensitivites will have to be improved by at least two 
orders of magnitude in future projects. This requires massive detectors operated by large inter- 
national collaborations. Accordingly, ASPERA has funded two projects, EURECA [9] based 
on solid state detectors, and DARWIN [10] using noble liquids. The latter consortium is to 
deliver the concept for a multiton LAr/LXe detector by 2013 which could become operational 
around 2016, reaching a final sensitivity of 10~^^ pb. Fig. [2] shows the possible layout of such 
a detector. 
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Figure 2: Left: Institutions participating in DARWIN. Right: Possible layout of a large under- 
ground detector using 8t ofLAr and 5t ofLXe, surrounded by water Cerenkov veto shields jjJJ^. 



2 Experimental sensitivities in noble liquids 

One often assumes that DM is distributed within a non-rotating spherical halo larger than the 
Milky Way. The velocity distribution is taken as Maxwell-Boltzmann with a dispersion of 220 
km/s. The solar system rotates around the galactic center with a velocity of 230 km/s and the 
earth around the sun with a velocity of 30 km/s. Taking into account the 60° inclination of the 
ecliptic one then finds that the average velocity of WIMPs entering a terrestrial detector is ~245 
km/s in June and ~215 km/s in December. The differential scattering cross section on a nucleus 
is given by fil2il 

da _ GlCF\q^) 

where is the momentum transfer to the nucleus, Gp the Fermi coupling constant, C a constant 
describing details of the interaction, F a nuclear form factor and v the WIMP velocity in the 
detector. The event rate in a detector can be calculated [13] by folding the cross section with 
the Boltzmann distribution, assuming a density of WIMPs of about 0.3 GeV/cm^. The cross 
section consists of two components, a nuclear spin dependent part and a spin independent part. 
The former is larger on single nucleons, but constructive interference enhances the latter by 
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a factor of for nuclei. Small WIMP velocities are favored by the cross section (eq. [B 
and the form factor suppresses large recoil energies. Hence the distribution falls off (roughly 
exponentially) as a function of nuclear recoil energy. This makes WIMP detection difficult due 
to the low energy background. Fig. |3] shows the recoil energy distribution in It of liquid argon 
and liquid xenon detectors assuming a cross section of 10^® pb. Recoil energy spectra in LAr 
and LXe are different. These liquids are therefore complementary in providing a crosscheck in 
case a WIMP signal is observed. With a nuclear recoil energy threshold of 30 keV, a WIMP- 
nucleon cross section of 10~® pb would yield 1 event per day and per ton of LAr. Even alt 
LAr detector could improve significantly the current upper limits, provided that the threshold 
of 30 keV and sufficient background rejection can be achieved. Statistics permitting, the annual 
variation could even be used as additional evidence for the detection of WIMPs. 
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Figure 3: Daily WIMP detection rate for a It detector as a function of detection threshold for 
LXe and LAr, assuming a WIMP mass of 100 GeV and a cross section of 10~^ pb. The top and 
bottom curves show the limits of annual variations due to the rotation of the Earth around the 
Sun (after ref ^T3\l). 



3 Light detection in a It liquid argon detector 

The ArDM Collaboration [fT4l uses the double-phase argon technique pioneered by the WARP 
collaboration [1151 . ArDM envisages a much larger and more ambitious detector (It of LAr). In 
contrast to WARP one plans to detect both charge and luminescence from WIMP interactions. 
Charged particles lead to ionization and excitation of argon atoms, forming excimers with the 
lowest singlet and the triplet excited states decaying by VUV photon emission in a narrow band 
around 128 nm. The singlet and the triplet states have different decay times, respectively ri ~ 5 
ns and T2 — 1.6 /us in liquid argon (LAr) [[T6l . The slow component is very sensitive to the 
purity of the argon Wl\ . Heavily ionizing particles such as a's or nuclear recoils contribute 
mostly to the fast (< 50 ns) decaying component, while the contribution of electrons and 7's 
to the slow component is larger. For minimum ionizing projectiles, such as electrons and 7's, 
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the component ratio CR (fast to total intensity) is ^ 0.3, while for a's and nuclear recoils 
CR ~ 0.8 [il6il . In addition, the ionization yield is much lower for nuclear than for minimum 
ionizing particles, due to quick recombination. Both features can be used in LAr to reduce 
background in WIMP searches. 

LAr has a density of 1.4 g/cm^, a refractive index of 1.24 and boils at 87 K. LAr generates 
about 5x10^ VUV photons/MeV for minimum ionizing particles (in the absence of electric 
field). LAr is cheap (1.5 CHF/i) compared to xenon, but has the disadvantage of containing the 
radioactive nuclide '^^Ar {(3 - emitter with a half life of 269 yrs and 565 keV endpoint energy). 
The activity of natural argon is 1 Bq/kg which induces a background rate of roughly 1 kHz in a 
It detector [fTSl . 




Figure 4: Left: Detector concept. Middle: PMT array and HV divider chain. Right: WLS-foils 
and PMTs under UV illumination / f20l/ . 

Details on the ArDM detector can be found in recent publications [T9l, l20l . Briefly, the 
working principle is as follows (Fig. |4l left): in LAr a WIMP collision leading to 30 keV nu- 
clear recoils produces about 400 VUV (128 nm) photons, together with a few free electrons. 
The latter are drifted in a strong vertical electric field and are detected in the gas phase by a 
large electron multiplier (LEM) above the surface of the liquid, while the VUV scintillation 
light is shifted into blue light by a wavelength shifter (WLS) and detected by cryogenic photo- 
multipliers at the bottom of the vessel. Fig. |4] (middle) shows a sketch of the detector with the 
400 kV HV divider. The detector and its recirculation and purification systems are currently 
installed in building 182 at CERN for preliminary readout tests. The detector for the charge 
extracted from the liquid into the gas phase with subsequent multiplication in a prototype LEM 
still needs to be installed. 

At the University of Zurich we have designed and built the light readout system for the It 
LAr detector. We measured the light yield dependence on the WLS thickness of evaporated 
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Figure 5: Left: disks of ESR-foils covered with TPB under ambient and UV light. Right: Light 
yield in photoelectrons (pe) as a function ofWLS thickness. 



specular reflectors foils. The measurements were made with a small cell containing gaseous 
argon at atmospheric pressure. Two materials were selected, Tetratex and ESR-foils from 3M, 
and the optimum thickness of the TPB layer on the reflector was determined experimentally 
ETI . The best conversion efficiency, uniformity and reproducibility were achieved by evapo- 
rating TPB on the reflector. Several disks of ESR foils (diameter 70 mm) were covered with 
TPB layers of different thicknesses. The TPB thickness was determined by weighing the disks 
before and after evaporation. The increase in brightness with TPB thickness is apparent under 
UV illumination (300 nm) (Fig. [51 left). The response to 128 nm light was determined using 
scintillation light from gaseous argon and an a-source. The light yield of the slow component 
(t2 = 3.2 /iS in gas ifTTl ) is shown in Fig. [5] (right) as a function of thickness. The data are 
consistent with a saturation of VUV conversion efficiency above 1 mg ■ cm^^. Accordingly, 
fifteen Tetratex WLS sheets (120 x 25 cm^) were coated with 1 mg ■ cm^^ TPB to cover the 
cylindrical side walls inside the electric field shapers (Fig. |4l middle). 

The light detection system consists of fourteen 8" hemispherical photomultipliers (PMT) in 
a staggered arrangement at the bottom of the vessel. We have investigated PMTs for their func- 
tionality and quantum efficiency at low temperature. The best result was obtained with Hama- 
matsu PMTs (R5912-MOD) manufactured with Pt-underlay. A light but sufficiently strong me- 
chanical support was constructed to withstand the buoyant force (^IkN) acting on the PMTs in 
liquid argon. The PMT glass was coated with a thin WLS layer of a transparent TPB-paraloid 
compound to increase the VUV light yield. Fig. |4] (right) shows the mounted foils and the 
PMTs under UV illumination. Many details can be found in a recent thesis [|22l . 

The ArDM detector was filled for the first time with It of LAr in 2009. Several important 
parameters such as stable cryogenic operation in high LAr purity, high scintillation light yield, 
and detection of events down to energies of tens of keV's could be verified. The test was 
performed with a partial light readout assembly consisting of half of the PMTs, no electric field 
and no charge readout. The LAr purity was found to be constant over the measurement time of 
three weeks by monitoring the decay time r2 of the slow component of the light signal [|20l . 

The measurements were done with external sources such as ^^Na, delivering positrons (an- 
nihilating into two 511 keV 7's) and monochromatic 1275 keV 7's. The light yield produced by 
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Figure 6: Light yield in the It LAr detector for 511 keV j's (in photoelectrons, p.e.) at two 
positions of the "^"^Na- source. The measurements are in black, the simulated data in red (from 
ref 12Q]). 



one of the 511 keV 7's, following (multiple) Compton scattering, was measured by triggering 
with a 4" Na(Tl) crystal on the second 511keV emitted in the opposite direction, and on the 
1275 keV 7. Fig. [6] shows the light yield distributions for the source located at two different 
vertical distances from the photomultiplier array (black bars) [|20l. The light yield (in photo- 
electrons, p.e.) was calculated by Monte Carlo simulation as a function of energy deposit, and 
the distributions (red bars) compared to the measured ones. Good agreement was found with 
an average light yield of typically 0.4 p.e./keV, which is roughly half of the yield that would be 
obtained with a completed detector (14 PMTs). Note that the light yield for nuclear recoils is 
lower due to quenching, typically 30% of that for electrons in the few 10 keV range, see Fig. 
[I2I below. We are therefore confident to be able to reach our goal of 30 keV threshold in It of 
LAr for WIMP detection. 



4 LAr response to neutrons 

The light yield of nuclear recoils in LAr are poorly known, especially below 50 keV (see e.g. 
ref. [23]). A good way to calibrate the light output and to study the response of LAr or LXe to 
nuclear recoils is to scatter a beam of monoenergetic neutrons of energy Tn on argon or xenon 
nuclei and to measure the light yield as a function of scattering angle 6, from which the recoil 
energy can be calculated according to the formula 



2T„ 



l + A- cos'^e - COS eVA"^ + - 1 



2TnA 



'1 — cos I 



(2) 



where A is the atomic mass (A ^ 1). The method is illustrated in fig. |7] (left) and the energy 
deposits are plotted in Fig. |7](right) as a function of scattering angle for T„ = 2.45 MeV incident 
neutrons. 

In fact, neutrons (mainly from cosmic ray induced spallation) are a serious source of back- 
ground in WIMP searches, due to the huge cross sections for neutron-nucleus scattering, many 
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Figure 7: Left: Principle of a scattering experiment with monoenergetic neutrons. Right: nu- 
clear recoil energies in argon and xenon for neutron energies of 2.45 MeV. 



orders of magnitude larger than the ones for WIMP- nucleus scattering. This background con- 
tribution can be reduced by operating underground and by shielding the detector with \ow-Z 
materials. However, the large neutron-nucleus cross section also means that the mean free path 
of neutrons is of the order of a few cm, leading to multiple scattering, quite in contrast to 
WIMPs which interact only once in the target. This gives an additional way to suppress neutron 
induced background, provided that multiple scattering can be detected from the enhanced light 
yield, and/or by identifying the multiple interaction vertices in the detector, e.g. with a charge 
detection device. 





Figure 8: Sketch of the neutron generator, the LAr cell, the cryocooler and the four liquid 
scintillators (LSCs) to detect scattered neutrons. The right picture shows the LAr cell with its 2 
PMTs on top and bottom of the LAr volume. 

We have therefore purchased a monoenergetic neutron source f24] based on the reaction 
dd — J-^He n, and are setting up a scattering experiment with coUimated 2.45 MeV neutrons. 
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The target is a small (<li) test cell (77 mm high and 74 mm in diameter) and liquid scintillation 
counters (LSC, EJ301 from SCIONIX) detect the scattered neutrons in coincidence as a function 
of scattering angle (Fig. [8]). To reduce the measurement time we use four LSCs to cover various 
angles in parallel. 

The fusion chamber is surrounded by a 90 cm diameter shield of borated polyester and the 
experiment confined within a radiation controlled fence in our laboratory at CERN (Fig. |9l left). 
Residual radiation (mainly from scattered neutrons and X-rays) is well below the authorized 
limit of 2.5 //Sv/h. The neutrons are coUimated through a polyethylene orifice within roughly 
1% X 47r sr. The neutron flux (up to 5x 10^ n/s in An, according to specifications) is controlled 
through the applied high voltage and discharge current. 




Figure 9: Left: 2.45 MeV neutron fusion generator inside its safety fence in our laboratory 
at CERN. The source itself is shown in the inset. Visible are also the X/^-ray monitor {top 
white box), the neutron radiation detector {black sphere) and the electronic rack {right) with 
the measuring electronics and the HV control units. Right: Neutron generator flux into 4% {see 
text). Straight lines are drawn to guide the eye. 

The neutron flux was measured with a 5" LSC located at the collimator exit. The polyethy- 
lene collimator was surrounded by a 2 mm thick lead box against X-rays, inserted into the 
polyester shielding. According to NSD-Fusion [24J, the highest possible voltage and current 
are 100 kV, resp. 15 mA, corresponding to a flux of 10''n/s into Air. The flux increases pro- 
portionally to the current and to the voltage V^^. Measurements were made with and without 
polyethylene collimator and the rates corrected for the solid angle, assuming a detection effi- 
ciency of the LSC of 80%. A NIM MPD-4 module was used to discriminate between neutrons 
and X-rays. Figure |9] (right) shows the measured neutron intensity (blue, red and orange dots) 
compared with the values specified by NSD-Fusion (green crosses for 10mA and orange crosses 
for 15mA [|24l). Insufficient lead shielding leads to a strong X-ray contamination. The measured 
flux (red dots) was roughly one order magnitude smaller than anticipated, with a maximum of 
10^ n/s. 

The neutron energy distribution (smeared by collimator scattering) must be known to mea- 
sure the light yield from LAr accurately. Fig. [10] (left) shows the neutron energy distribution 
measured by a 5" LSC placed at the exit of the collimator. Ideally the spectrum should be flat 
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(without collimator scattering, with infinite resolution and single n-scattering). The spectrum 
of Fig. [10] (right) was obtained by unfolding the response of the LSC to neutrons from a ra- 
dioactive source. The response was obtained with an AmBe-source, by measuring the neutron 
energy through the time-of-flight between the source and the neutron counter over a distance 
of 1 m. The start time was determined from the 4.4 MeV 7 detected in a BGO crystal located 
close to the AmBe-source. 




Integrated pulse hight [nVs] Energy [MeV] 



Figure 10: Left: Pulse height distribution induced by neutrons in the LSC at the exit of the 
collimator. Right: Incident neutron energy spectrum obtained after unfolding the response of 
the LSC. 




We have started to measure the scintillation response of LAr to nuclear recoils with the 
LAr cell shown in Fig. [8] (right). Wavelength shifting reflectors (Tetraphenyl-Butadiene, TPB, 
on Tyvek foils) were mounted on the inner walls of the cell to convert the 128 nm light into 
400 nm. The cell was read out by two Hamamatsu R6091-01MOD PMTs with Pt-underlay 
for cryogenic operation. An internal ^^°Pb-source emitted 5.3 MeV a's and up to 1.2 MeV 
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electrons. We recall that the component ratio Ci? is the ratio of integrated light yield during 
the first 50 ns to the total light yield. Thus a high CR corresponds to the emission of light 
with mainly the fast component. As mentioned above, heavily ionizing particles such as a's or 
nuclear recoils lead to a large CR value. Fig. [TT] shows the component ratio CR from one of 
our first measurements of argon luminescence with the fusion generator. A clear contribution 
from neutron induced nuclear recoils is observed at CR ~ 0.8. 

A first measurement of the scintillation efficiency for nuclear recoils relative to electrons was 
performed with a 5" LSC at 65° from the incident beam direction, at a distance of 50 cm from 
the LAr cell. The time-of-flight between the LAr cell and the LSC could be determined off-line 
and used to remove background, e.g. from multiple neutron scattering in the cell. The reference 
time was determined with a Na-source located at equal distance from the LAr-cell and the LSC, 
and using the two back-to-back 5 1 1 keV 7's. At 65° the argon recoil energy is 69 keV with 2.45 
MeV incident neutrons. This corresponds to a time-of-flight of 23 ns for 2.38 MeV neutrons 
flying to the LSC. Good pulse shape separation between proton and background (7-induced 
electron recoils) could be achieved with the LSC by analogue pulse shape discrimination. 




Nuclear Recoil Energy (kev ) Integrated pulse height [nVs] 



Figure 12: Left: Measurement of the relative scintillation efficiency in LAr as a function of 
recoil energy (from ref. / f23l/ ]. Our preliminary measurement is shown by the red data point. 
Right: Energy spectrum from Am-decay measured with the refurbished LAr cell, showing the 
60 keV line. 

Several cuts were applied in the offline analysis. For example, we required the component 
ratio CR to be larger than 0.6 (see Fig. \TT}. We also rejected events close to the PMT windows 
by comparing the signals from the two PMTs. A time-of-flight window between 18 and 31 ns 
was selected. We obtained a scintillation efficiency of 0.27 (with about 25% error) in this first 
attempt. Fig. [12] (left) shows our data point compared with data from ref. [|23l . 

Light yields of nuclear recoils are usually determined relative to electronic recoils. The 
electronic light yield is determined with various external 7-sources and also with a ^^Kf'" source 
which can be connected directly to the gaseous phase in our setup. 

To improve the light yield and the measurements with radioactive sources we upgraded our 
LAr cell to reduce the thickness of the stainless steel vessel and the LAr volume. Tetratex foils 
(from Donaldson Membranes) were used instead of the Tyvek reflectors, and TPB-Paraloid 
coating was replaced by 0.08 mg-cm^^ of evaporated TPB. Fig. [I2](right) shows the spectrum 
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of the 60 keV line from Am-decay. During data taking the mean life of the slow component 
was much lower than the established value of 1.6 fis, due to the impurity of the LAr (see e.g. 
ref. [|20l '). By extrapolating the light yield to maximum purity we could set a lower limit of 3.2 
p.e./keV. 

The setup has now been improved by the addition of a cryocooler (Gifford-McMahon-type) 
mounted on top of the LAr cell. A cooling serpentine is bonded to the temperature regulated 
cold head and provides the liquefaction of the recirculated gas. Gas purification is achieved by 
two cleaning cartridges. They reduce the O2 and H2O contamination below 20 ppm. A further 
upgrade of the LAr cell with a smaller active volume (and less surrounding material) to reduce 
multiple scattering and enhance the light yield is in progress. 

5 Conclusions 

The current experimental sensitivity on WIMPs is around 10^*^ pb for masses of ~100 GeV. 
The goal of future experiments using solid state detectors or noble liquids such as LAr or LXe 
is to reach 10^^° - 10^^^ pb, a sensitivity within theoretical predictions from supersymmetry. 
The goal of the DARWIN consortium is to deliver the design for a LXe/LAr detector to become 
operational around 2016. Liquid argon technology is still in its infancy and intensive R&D is 
ongoing, in particular to determine the light yield from heavily ionizing nuclear recoils. For the 
first time a large It LAr argon prototype detector (ArDM) was operated on the surface with a 
light yield of 0.8 photoelectrons/keV electron equivalent. These are encouraging prospects to 
detect WIMPs in LAr with a threshold below 50 keV. 
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